arXiv:nucl-th/9706039vl 17 Jun 1997 


Antisymmetrization in the Multicluster Dynamic Model of Nuclei 

and the Nucleon Exchange Effects 

R.A. Eramzhyan^’^, G.G. Ryzhikh^’^, and Yu.M. Tchuvirsky"^ 

^ Institute for Nuclear Research of the Russian Academy of Sciences, Moscow 117312, Russia 
^ Gyclotron Institute, Texas A&M University, Gollege Station, TX 77843-3366,USA 
^ Faculty of Science, Northern Territory University, Darwin NT 0909, Australia 
Institute of Nuclear Physics of Moscow State University, Moscow 119899, Russia 


Abstract 

A modified version of the Multicluster Dynamic Model of nuclei is proposed to 
construct completely antisymmetrized wave functions of multicluster systems. An 
overlap kernel operator is introduced to renormalize the total wave function after 
antisymmetrization between nucleons in different clusters. A group-theoretical 
method is developed to analyze the role of the exchange effects arising in the 
calculation of the various observables of multicluster systems due to this antisym¬ 
metrization. 

The Antisymmetrized version of the Multicluster Dynamic Model is applied 
to the six-nucleon systems treating them as a-2N ones. The static and dynamic 
characteristics of the six-nucleon systems manifested in electron and vr-meson 
scattering, muon capture, /3-decay, pion photoproduction, etc., are calculated. 

Significant progress is achieved in describing of variety of dynamic observables 
of the six-nucleon systems as compared to the multicluster dynamic model. In 
most cases calculated static and especially dynamic characteristics are in a good 
agreement with the experimental data. 

Keywords: nuclear structure, nuclear models, light nuclei, nuclear reactions, electron scatter¬ 
ing, pion scattering, pion photoproduction, muon capture. 


1 


PACS numbers: 21.60.Gx, 27.20.+n, 25.30.-c, 25.80.-e 


Typeset using REVTgX 


2 



I. INTRODUCTION 


To describe the structure of light nuclei (A>4) in great detail very sophisticated approaches 
are now being elaborated. Recently, Green’s Function Monte Carlo calculations have been re¬ 
alized for nuclei with A< 6 starting from the realistic nucleon-nucleon interaction. At the 
same time the stochastic variational method (this method was hrst applied in calculations of 
nuclear structure within the multicluster dynamic model (MDM) in Ref. [Q]) has been extended 
to the systems with A< 10 treating all nucleons variables on the same ground (see, for example. 
Ref. 1^). Before that the K-harmonic method and the extended shell model have been 
used to learn more about the structure of the Ip-shell nuclei. However, practical applications 
of all these methods is not an easy task. That is why there are not many observables analyzed 
with these newly developed methods. For this reason they are unable to displace the traditional 
approaches, considering the fact that the latter are improved as well. 

The Microscopic Cluster Model (MCM) still continues to play a role as a very successful 
approach well suited for description of the structure of light nuclei. This model in its develop¬ 
ment has passed many stages, starting from its simplest realization within the framework of the 
single channel Resonating Group Method (RGM) P| and coming hnally to the many-channel 
multicluster MCM Mfl. 

This model proceeds from the tendency of the nucleons to form clusters inside the light nuclei 
due to their strong binding. In the MCM the A-nucleon wave function is approximated by a 
superposition of the antisymmetrized basis functions in which nucleons have been grouped into 
some clusters — a, ^He, ^H, d etc. The cluster internal wave functions are supposed to be known 
from the beginning but the relative motion function is obtained from the solution of the dynamic 
equation (or also from the system of coupled equations) derived straightforwardly from the initial 
A-nucleon problem. 

This approach in its hnal version enables us to incorporate into the calculation various mul¬ 
ticluster conhgurations, to take into account some effects of the cluster excitations and their 
rearrangement, to treat nucleons on the same ground as clusters, etc. Finally a consistent and 
successful description of bound, resonant and scattering states of many nuclei P-[T^ has been 


3 





achieved. 


At the same time extensive applications of the MCM have revealed some dehciencies which 
show up mainly in the case when the short-range correlations among nucleons appear to be im¬ 
portant. Rehned realization of the model requires that more realistic nucleon-nucleon interaction 
should be used and, consistently, more realistic wave functions of clusters should be involved in 
the calculation. But such a program is extremely difficult to realize in practice. Indeed the 
cluster-cluster potential which finally arises in the MCM comes about from the folding of the NN 
potential over the intrinsic wave functions of the constituent clusters. This potential appears to 
be very complicated, highly nonlocal and energy dependent due to the exchange effects and the 
composite nature of clusters. If one starts from the cluster internal wave functions with both 
short range and tensor NN correlations and uses the realistic NN forces with the noncentral com¬ 
ponents, then one arrives to such a complicated cluster-cluster potential that it becomes almost 
useless for practical applications. That is why a simple intrinsic wave function of cluster and a 
simplified effective NN forces are usually used in practice. Of course in this case this approach 
looses its ’’pure” microscopic character and it becomes model dependent. 

On the other hand, when one starts from the NN potential with a strong short-range repulsion 
and uses noncorrelated wave functions of clusters then the results of the calculation become highly 
dependent on the number and the type of the cluster channels taken into account In 

spite of this fact, many calculations are still undertaken in this approximation with forces like 


the Hasegawa-Nagata or the Eikemeier-Hackenbroich one (see, for example. Refs. |1l5|,pT 


There is another problem which is a principal one for the MCM. Derived from the more or less 
realistic NN potential, the cluster-cluster potential of the MCM does not reproduce accurately 
the phase shifts for free clusters. This problem was studied thoroughly in the case of the simplest 
cluster system, namely for the a-N (see Ref. [|^). Here the N-a potential was constructed 
starting from a realistic NN potential and a very complicated a-particle wave function. It turned 
out that even in this simple case, in order to reproduce the experimental a-N phase-shifts the 
starting realistic nucleon-nucleon potential had to be modified significantly. It means that in 
the MCM, a phenomenological (effective) NN-potentials should be used rather than the realistic 


ones. 
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At the same time a microscopically inspired Potential Cluster Model, which is also called the 
Multicluster Dynamic Model — MDM, has been also applied successfully for the calculation of 
the structure of light nuclei. Formulation of the MDM is given, for example, in Ref. [|^]. In this 
approach the relative motion function of clusters is obtained as a solution of the fc-cluster dynamic 
equations with the cluster-cluster potentials derived directly from the phase-shifts analysis. Thus, 
in the MDM the on-shell behavior of the interaction amplitudes are taken from the outset with 
higher precision as compared to the MCM. 

In the MDM there is no problem at all with the intrinsic cluster wave functions due to the 
factorization of the matrix elements squared for the nuclear observables. Due to this factorization, 
it is not the intrinsic cluster wave functions but rather the magnitudes of their observables (that 
are directly taken from experiments) which are used here. In this way many observables of 
the clusterized nuclei were obtained in a good agreement with the experimental data (see, for 


example. Refs. ||T9|-^). Another attractive point of the MDM is its much more simple formalism 
as compared to the MCM. These simplihcations come from the approximate treatment of the 
Pauli-principle for the nucleons from different clusters. 

However, there are some observables where just the nucleon exchanges between different 
clusters are of primarily importance. That is why we have modified the MDM to construct dually 
completely antisymmetrized wave functions. This modified version of the model was called the 


AMDM — the Antisymmetrized version of the MDM. It was shown in Ref. that by using the 
AMDM we extend the list of well reproduced observables as compared to the MDM. In particular, 
it was possible to describe all longitudinal and transversal elastic and inelastic electromagnetic 
form factors of ®Li simultaneously for a broad interval of the momentum transfer. At the same 
time the AMDM formalism is still more simple and transparent than the MCM one. 

The method of construction of the totally antisymmetric wave function from the MDM func¬ 


tion which had been used in our earlier works p4p3| has some limitations. The fact is that after 
antisymmetrization of the MDM functions between nucleons from the different clusters, the set 
of total wave functions obtained loses the orthonormalization property. First of all one has to 


restore the total normalization of the wave functions. In Ref. done by the simplest way 

of multiplying the wave function by a constant. This is a common procedure for the shell model 
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approach where the number of oscillator quanta for each wave function is hxed. However, in gen¬ 
eral case such a procedure is not quite correct. The reason is that the constant renormalization 
(C-renormalization as we will abbreviate it further) does not restore the orthogonality between 
different antisymmetrized wave functions with the same total quantum numbers. Moreover, this 
C-renormalization changes the asymptotic part of the relative motion wave function where the 
antisymmetrization must have no influence at all. 

Experience accumulated in the MCM suggests that the renormalization with the integral ex¬ 
change kernel (it'-renormalization), which appears primarily in the RGM, will allow us to solve 
this problem in the general case as well. It is very important that there is a kernel renormal¬ 
ization that allows to determine the class of observables which are exactly conserved under the 
antisymmetrization [^. Some aspects of the renormalization problem have been discussed al¬ 
ready in Ref. for ®Li which was treated as a three-cluster system. In this paper we present 
new precisely renormalized version of the AMDM and its application to many observables of the 
six-nucleon systems. The six-nucleon system is particularly singled out in nuclear physics and 
has been studied in many papers as a three-nucleon system. Six nucleons in nuclei exhibit quite 
properly the typical properties of all light nuclei, such as the decisive role of the Pauli principle, 
the substantial spin-orbit splitting of levels, a signihcant greater role of the interaction of valent 
nucleons with core compared with the interaction among them, the important contribution of the 
short-range correlations of nucleons, etc. These circumstances, as well as numerous experimental 
data obtained for six-nucleon systems including the strong, electromagnetic and weak interaction 
processes permit their theoretical studies to be regarded as a theoretical laboratory for nuclear 
physics in a broad sense allowing to test new theoretical approaches. 

This paper is organized in the following way. In Section II the formalism of the AMDM is 
briefly described and a group-theoretical analysis of matrix elements for various observables is 
performed. In Section III many observables of the A=6 nuclei treated as a three-cluster a-2N 
systems are analyzed within the framework of both versions of the AMDM — with constant 
and kernel renormalization. All these results are compared to those obtained within the MDM. 
Besides the traditional observables — static properties of ground state, beta-decay and electron 
scattering — we are discussing muon capture, pion scattering and pion photoproduction. All 
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these reactions proceed via similar nuclear transition operators. Studying all of these transitions 
does not simply lead to a large amount of redundant information, but, rather, it generally leads to 
complementary information. So the most informative approach for studying the nuclear structure 
is to carry out a simultaneous study of all of them. We conclude in Section IV. 

II. GENERAL PROPERTIES OF THE AMDM 

A. Construction of the MDM wave function 

The AMDM is essentially based on the MDM. Therefore it makes sense to run briefly through 
the basic ingredients of the MDM. 

Let A identical fermions (nucleons in our case, but they could be quarks, etc.) be distributed 
over k < A clusters. Individual nonclusterized (valent) fermions, if they appear in this case, are 
treated as clusters as well. 

In the MDM, for a system made up of the k clusters the wave function, with given total 
momentum J, total isospin T and their projections M and is dehned as 

({?"})= i: (') 

where Jc is the channel spin, denotes the complete set of the intercluster Jacoby coordinates. 

'Llmz,({'C}) is the function of the cluster-cluster relative motion with the total angular momentum 
L. It depends on the corresponding set of the Jacoby coordinates Each coordinate is 
divided into angular and spatial parts: 

We hrst dehne some notations to be used in this Section. Let’s x to be any momentum (/ — 
angular momentum, s — spin, j = I + s — total momentum, t —isospin) and rrix its projection, 
then the symbol will denote the coupling scheme of the individual momenta Xi and their 

projections to the intermediate Xi2) ■■■Xi2...k momenta and the total, x, momentum and its 
projection: 

Correspondingly, the symbol with double down index is dehned as: 


7 



l-.V = l-.V l-.V . 

I Jxmx,ymy I j xrrix I j yruy 

The following abbreviations are used to list the intermediate momenta in eq. (0): 

{j} = {jl2, jl23, ■■■, jl2...k, Jc), 

{t} = iti2,ti23, 

{/} = (/i, I 2 , h-l, h2, ^1235 ^12...A:-l5 L), 


( 3 ) 


(4) 


where k is the orbital momentum conjugated with the Tth Jacoby coordinate of the relative 
motion; L is the total angular momentum of the system, Jc is the spin of the channel. The total 
momentum J equals to the sum J = Jc + L. Finally the complete set of quantum numbers {a;} 
which characterizes any channel of the system is defined as 

= (5) 


According to the MDM, the internal wave function of the clusters is not modihed in the 
nucleus and it is built up from the wave functions of free clusters in their ground states: 

k 


^int ^ JJ- ^ 




,({C'}). 


( 6 ) 


2=1 

where is the set of the internal Jacoby coordinates for the i-th cluster and ji and ti 

are their total momentum and isospin with projections m* and mt-, respectively. The role of 
the cluster polarization inside a nucleus has been studied recently within the MCM. It was 
shown in Refs. ||T3| , p!^ that the a-cluster is polarized very weakly inside light nuclei. Therefore, 
the assumption that the a-cluster can be considered as a free a-particle is rather good in any 
multicluster approach. 

The wave function of the cluster-cluster relative motion is decomposed into the 

angular F{z}m,( 0 fhe radial parts. The angular part of the relative wave function 

can be written as: 


'k-l 


noAf,(o= n Ki 


. 2=1 


k-1 


LMl 


(7) 


It is very convenient to present the radial wave function of the relative motion in terms of 


the superposition of gaussians (it was demonstrated in many papers, see for example Ref. |T^) 
as follows: 
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( 8 ) 


fe-i 


= ECJ“> n l6l'*exp(-4"'??) , 


2=1 


where and are the expansion coefficients. The set of these coefficients together with 
the set of the ’’channel” quantum numbers {a;} determine completely the MDM wave function. 
The variational method for solution of the /c-cluster problem is well suited for hnding these 
coefficients. As it was shown in Refs. [^,^], such an approach allows to obtain a very accurate 
solution of the nuclear few-body problem. It is the gaussian parameterization of the MDM wave 
function that makes it possible to get analytical expression for most matrix elements both in the 
MDM and in the AMDM (see, for example, Refs, ). 


As it was shown by Saito that sophisticated nonlocal energy-dependent ’’microscopic” 
cluster-cluster forces appearing in the RGM can be successfully approximated by a simple local 
A'-independent ones if they are supplemented by the orthogonality requirement based on the 
Pauli principle. In such an approach, both the phase shifts and the nodal structure of the RGM 
wave functions can be quite well described in a very simple way. This is the basic approximation 
for the Orthogonality Gondition Model (OGM) pTjpSl and, in fact, for any MDM. In the version 
of the MDM to be used throughout this paper, the Pauli-forbidden components are excluded 
from the solution of the MDM dynamic equations by means of the pseudopotential technique 
1^ , 1^ . As we already stated before, such ’’microscopically inspired” intercluster forces have 


even some advantage compared to the ’’true microscopic” ones. This advantage consists in 
more accurate description of the low-energy cluster-cluster phase shifts and, therefore, the most 
important on-shell properties of the cluster-cluster forces. 

However, exclusion solely of the forbidden components does not account for all the conse¬ 
quences of the Pauli principle for the multicluster system. Exchanges between nucleons in the 
different clusters are completely omitted in the MDM. That is why some types of observables 
can not be explained in the framework of this model. Therefore we came to recognize that the 


totally Antisymmetrized MDM (the AMDM) should be used in general case Notice that 

the AMDM is, in fact, a modihcation and generalization of the OGM. 
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B. Wave functions of the AMDM 


In the initial version of the AMDM [Q named here as AMDMc the total wave function was 
constructed by means of direct antisymmetrization of the MDM wave function: 

1 


4 .. = 


MDM 


(9) 


As we have noted above, in this case we use a constant Q to normalize the total function 
That is why we call this version of the model as the AMDMc — Antisymmetrized Multicluster 
Dynamic Model with the renormalization to the constant. A in eq.(^ is an antisymmetrizer: 


i = (^1 + H-irpj , 

where the normalizing factor D has the following form: 

A! 


( 10 ) 


n = 




( 11 ) 


Here Aj is the mass of the individual cluster i, P is the operator permutating nucleons between 
different clusters and p is the parity of such a permutation. 

In a new version of the AMDM, named the AMDMk, the overlap kernel operator K is 
introduced to renormalize the total wave function after antisymmetrization: 






( 12 ) 


This operator K acts in the Hilbert space of the (3/c — 3)-dimensional orbital functions of the 
/c-cluster system (cluster subspace). It acts on the function in the following way: 


«({?})= /<l{f}A'({f'.f })*(«'}). (13) 

The overlap kernel is dehned as a projection of the antisymmetrizer onto the 

cluster subspace of A-nucleon variables: 

}) = a . n '^(^i - I -41 n '*(& - r*)) ■ (11) 

\ ’ i=l ’ i=l / 

The eigenfunctions aiid the eigenvalues of this kernel satisfy the usual equation: 
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( 15 ) 


and 




( 16 ) 


for 


^ 0. (17) 

In the MDM, any nnclens is nsnally considered as a system composed of some nnmber of 
magic clnsters (a-particles, etc.) in their gronnd states and n-valent nncleons. In this 
Section, to make the gronp-theoretical analysis more transparent we restrict the general case by 
the following two conditions: 

1. The magic clnsters are consider to be the SU(3) and SU(4) scalars. So, in their wave 
fnnctions we neglect all components violating these symmetries. If the system is bnilt np 
from different clnsters (a-particles and for example) we will neglect also the difference 
in their oscillator parameters huj in order to conserve the SU(3) scalar character of 

2. The valent gronp of n nncleons is considered either as one light cluster with three or less 
nucleons (with the common oscillator parameter a;) or as nonclusterized nucleonic system 
if n < 2. 


These conditions lead to that the total spin S = Jc and isospin T are uniquely determined by 
the Young scheme [f] and Yamanouchi symbol (r). These quantum numbers are precisely those 
that characterize the MDM functions. 

As a consequence of the SU(3)-symmetry, the eigenfunctions of the operator K are the many- 
body oscillator functions with a dehnite number of oscillator quanta N and a hxed value of the 
Elliott symbols (A/r). Since the same Elliott symbol can characterize several eigenfunctions, an 
additional quantum number i is introduced to distinguish among them. The eigenfunctions of 
the operator K are characterized by the total angular momentum L. It is interesting to note 


that its eigenvalues do not depend on L at all. Finally for a complete definition of eq. ([l5D we 
specify a symbol v as: 
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p = {S,T,L,N,{\n),i). 


( 18 ) 


According to these restrictions the internal part of the MDM function can be written down in 
the factorized form: 



k 

Jc A/c 1 Mt 



(19) 


where <hj„t is a scalar spatial function and the rest is the spin-isospin one. To simplify the reading 
of the subsequent equations we will use the following notations for the spin-isospin function: 

( k 


XST - X[f]{r) - S n A 


Sjti 


( 20 ) 


2 = 1 


SMs,TMt 


Before concluding this subsection we note that these restrictions are used only in Section 
II for the sake of the formal analysis of the antisymmetrization effects. As to the numerical 
calculations in Section II, more complicated wave functions were actually used. On the other 
hand, it appeared that the results of the formal analysis are valid even in the case when the 
conditions 1 and 2 are not fulhlled exactly. 


C. Formal analysis of the matrix elements within the framework of the AMDM 

Before starting any discussion of the results of the calculations within the framework of the 
AMDM it is very instructive to make some classihcation of the nuclear matrix elements. For this 
reason let us expand the relative motion functions of the MDM onto the eigenfunctions of 
the overlap kernel: 

't({f})=i:w{f})- (21) 

Generally, the wave function in the multicluster approach can contain some admixture of 
the Pauli-forbidden components with = 0. After antisymmetrization these components have 
to disappear. To simplify the subsequent expressions, it is assumed below that such forbidden 
components have been eliminated from the beginning (i.e. Q ii = Q). So one can write 

down that: 

($ I 6 I $0 = E 

1 / 11/2 
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( 22 ) 


(i{*i„i*^,({fi})xi/iw} I O I -4{*mi4«({el})xj/.](,T}). 

For any operator O which is symmetric over the nucleon permutations one can replace an 
antisymmetrizer, say in the left hand side of eq. by the factor Let us then add to this 
expression the sum over the complete set of the nonantisymmetrized wave functions: 

1 = E I 'S‘inh.({&})xiri(.-))(4’.‘«h-({?l})X[ri(.-)>. (23) 

0>[/”]h”) 

where the sum over (j) includes all orbital excitations of the constituent clusters. Because the 
antisymmetrizer does not change the values of the total spin S, isospin T and Young scheme [f] 
of the system one arrives to the following expression for the matrix elements: 

($ I o I $') = f] ^ I O I 

i2iV2u{j)u(r”) 

I A I 4S4„,({?1})xot(,.|). (24) 

Here the index (j) = (0) labels the intrinsic orbital functions of clusters in their ground states. 
From eq. (p^ ) it follows that if the operator O is diagonal in the space of the orbital functions 
}) over the Yamanouchi symbols (r) then the matrix elements of the MDM and 
AMDM are equal to each other because in this case (j) = (0) only and according to expressions 
(0), (0 and (0) the matrix element of the antisymmetrizer reduces to Such a situation 

occurs when: 

1. An operator is of pure spin-isospin type. 

2. An operator is of particular spatial type being a function of the Casimir operators of the 
SU(3) group and its subgroups. 

Indeed in the last case this operator is diagonal because the wave functions in (P^ have a definite 
SU(3) symmetry, and the Elliott symbol (A/i) is a good quantum number for the system. 

It is interesting to analyze the difference between the matrix elements of the full tensor 
operators and the matrix elements of the pure spatial ones. For this purpose let us use the 
following representation of the wave function: 
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( 25 ) 


A I 'I>!S4^({S})X|/|(r)) = -ECi? I 

Uf ^ 

where ClP is the Young operator which projects the wave function onto the wave function with 
the dehnite Young scheme [/] and Yamanouchi symbol (r) (rij is a dimension of the corresponding 
representation): 

( 26 ) 

Here Dl{](P) is the matrix of the corresponding representation of the permutation group. 

After substitution of eq. (^5l) into eq. (H), one arrives at the following expression: 

(K I o I «>') = E KK ({?l})xi/l(r.) I 6 I 4S4>4{?i})\'i/']M) 

oy«y{<i>™ii>.-({S}) I df I (27) 

where 

Q. = (hS4.({6}) I df I dS4>.({6}))5 = S2-‘(>VEjt (28) 

Qy is the normalization factor of the wave function CH'^ \ 

Now we are ready to establish some relation between the renormalization factor Q in eq.@, 
which have appeared in the AMDMc and the eigenvalues of the overlap kernel 6^,: 

d = ^E(Md. (29) 

ii zy 

Let us compare now matrix elements of a pure spatial operator with matrix elements of a full 
tensor one. It follows from eq. that in the latter case only the terms with the Yamanouchi 
symbols (r) ^ (ro) contribute to the AMDM matrix element, contrary to the case of a pure 
spatial operator. Evidently such terms are not present in the MDM case independently of the 
tensor structure of any operator. That is why the matrix elements of a full tensor operator are 
expected to be modihed much more than the matrix elements of a pure spatial one when the 
exchange terms are taken into account. 

Eq. (1^) allows to establish some properties of matrix elements appearing in various reactions 
as well. To demonstrate this let us consider, for example, the cluster knock out reaction at large 
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momentum transfer, though it is beyond the scope of this paper. Nevertheless, it is instructive 
to say some words about this case. The amplitude of the corresponding process is proportional 
to the matrix element given by the expression (EZD- As follows from this expression not only the 
(j) = ( 0 ) term, but also the terms with (j) 7 ^ ( 0 ) contribute to the amplitude of this reaction. 
This fact was known in the theory of cluster knock out reactions for the special case of oscillator 
shell model wave functions. These (j) 7 ^ (0) terms correspond to cluster deexcitations according 
to Ref. B- Expression (^7l) demonstrates that the cluster deexcitation terms appear in the 
amplitudes of the cluster knock out reactions and similar processes in the general case as well. 
Despite the fact that in the case of cluster knock out one deals with the hnal state wave function 
which was named as an asymptotic one (because in the corresponding wave function components 
with the eigenvalues ~ 1 dominate in the expansion (^T]) ) the sum over (j) in (|^) is not 
reduced to the trivial term with (j) = (0). To demonstrate this let us compare two expressions. 


In the hrst one the antisymmetrizer is removed from the right hand side of eq. (|^) , in the second 
one — from the left hand side. It is easy to show that they will be equal to each other if the 
following relation holds 


Qui Qu2 ^^' 2 ) ( 30 ) 

for all terms with 6 ,^ 7 ^ 0. However, this is not the case due to the contribution from the terms 
with 7 ^ 1. So, the transition amplitudes from the bound cluster state, even to the asymptotic 
one, contain the deexcitation terms. In most calculations these terms were not taken into account 
at all. Some examples of the importance of taking into account this type of exchange terms are 
given in Refs. 0,0]. 

Finally, as follows from our analysis, all operators can be subdivided into three groups ac¬ 
cording to the magnitude of the exchange effects in their matrix elements: 

1. Full tensor operators. Their matrix elements are affected maximally when going from the 
MDM to the AMDMk. 

2. Pure spatial operators. Their matrix elements are less sensitive to the exchange effects as 
compared to the hrst case. If the Pauli forbidden components in the MDM function are 
eliminated in a thorough way then the exchange effects usually become rather small. 
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3. Operators whose matrix elements strictly conserved when going from the MDM to the 
AMDMk- They are either pure spin-isospin operators or pure spatial operators being some 
functions of Casimir operators of SU(3) group. 

In most cases the nuclear observables are associated not with a single tensor operator but with 
their sum. For example, an operator of the magnetic moment has two terms. One comes from 
the magnetization current and the other from the convective one. If the total angular momentum 
of the system equals to zero then the contribution of the last term vanishes. More generally, if 
the spin and isospin of the system under discussion are equal to zero both in the initial and in the 
hnal states, then any operator of physical observable appears as an effective spatial operator. In 
such a case the nucleon exchange effects should not be too large. On the other hand operators of 
some nuclear observables (for example the operator of Ml transition) at low momentum transfer 
can be considered with a good approximation as an effective spin-isospin one. In this limit the 
exchange effects for such matrix elements should be very small. 

By this discussion we hnish the group-theoretical analysis of nuclear matrix elements. In the 
next section we will discuss the results of the numerical calculations of the various observables 
of six-nucleonic systems treating them as a-2N systems. 


III. SIX-NUCLEON SYSTEM WITHIN THE FRAMEWORK OF THE AMDM 


In this section we discuss the AMDM application to the six-nucleon systems (®Li and ®IIe) 
treating them as the a-cluster in its ground state with L = S = T = 0 and two outer nucleons. 
The Jacoby coordinates for this system together with the conjugated angular momenta are shown 
in Fig.l. 

As an input to the AMDM calculations, the three-body MDM functions T (p, r) named as 


FUNCTION-92 in Ref. are used. These wave functions were calculated with a large vari¬ 
ational basis which includes all important (i.e. with the weight P>0.1%) components with the 
partial angular momenta /, A < 4. 

As has been stated in the Introduction, not the intrinsic cluster wave functions but the 
corresponding experimental data for the studied observables are used in the calculations within 
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the MDM. However, when going from the MDM to the totally antisymmetric wave functions one 
has to use explicitly the intrinsic wave functions of clusters in analogy with the MCM. 

The major part of the calculations within the MCM was performed using the simple OS"^ 
harmonic oscillator shell model function for the a-particle. Of all possible nucleon correlations 
only the breathing mode of the a-particle or the ^H-p (^He-n) channel was incorporated into the 
most advanced calculations (see, for example. Ref. iTl| ,[I^ , |3^ ). For many static observables, such 
a simplihed version of the a-particle wave function does not affect the hnal result. However, 
when large momentum is transferred to a multicluster system, then the short-range correlations 
inside the a- particle become very important. 

Unfortunately, it is known that the a-particle wave function calculated within a realistic dy¬ 
namic approach does not reproduce its Coulomb form factor at intermediate momentum transfer, 
when evaluated in the impulse approximation. To reproduce the experimental data either a very 
large contribution of the isoscalar exchange current has to be assumed or the existing wave 
function should be modihed to reflect the more complicated dynamic of this system. It seems 
that such a large contribution of the MFC in the a -particle, and therefore in the core of the 
six-nucleon system, is physically not adequate. This is the reason why we decided to modify 
the a-particle wave function and did it in a purely phenomenological way. When constructing 
this wave function we kept in mind that it should have the form convenient for performing the 
calculation of matrix elements in six-nucleon system with the antisymmetrized wave function of 
eq. (0). 

Following these requirements a phenomenological a-particle wave function was constructed 
in Ref. |2^ : 


* = exp ( -^ i:(rj - + C exp I i:(r, - Sf - ^(f. - S) 


N 


ai 


a2 


(31) 


k=l 


i=l 


k^i 


with the parameters ai = 0.6144/m“^, a 2 = 6.967/m“^ and C = —0.4506. Within the impulse 
approximation, this function hts the a-particle charge form factor in a wide region of momentum 
transfer. This function was used in Ref. 12211 to calculate the ®Li form factors within the AMDM. 


A signihcant improvement of the theoretical description of ®Li form factors at medium and high 
momentum transfer was achieved by using this modihed a-particle wave function. 
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The other advantage of the function (|^) consists in its similarity to the oscillator OS^ one. 
Their overlap (< T|Tosc >) is larger than 0.9. Therefore one would expect that the qualitative 
results of the group-theoretical analysis discussed above will be valid in this case as well. 

So the model with antisymmetrization is completely dehned and we are ready to start pre¬ 
senting the results of calculations within the model. We will discuss below the following three 
versions of the model: 


— the MDM, 

— the AMDMc , where the normalization with a constant is adopted, and 

— the AMDMk , where the integral kernel K is used for the renormalization. 


A. Operators and matrix elements involved in calculation of the observables of the 

six-nucleon systems. 

Three types of tensor operators are involved in our calculations of the various observables of 
the six-nucleon system. They are: 

^ TttJiqrj) '^{knikwm^ : umu)akmu{j)ywmSrj), (32) 

= intJiqrj) : umu)Y^m^{rj)pim{j) (33) 

and 

= *OoT“’“^(r,)(a,p;). (34) 

Here p is the nucleon momentum inside the nucleus and q is the momentum transfer to the nucleus; 
To and (To denote the unit operator (A: = 0,t = 0) in the isospin and spin space, respectively, and 
Tiz and aiz are the ordinary isospin and spin operators {k = l,t = 1), respectively; f{qrj) is the 
observable scalar function associated with the process which is being considered. 

We designate the nuclear reduced matrix elements associated with the above given operators 
in the following way: 
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< >= Mat , 


(35) 


< J,T,\\OZ-“-{pjmT, >= |lu.t.;p]AT 


(36) 


and 


< >= [Owu;p] 


AT 


(37) 


They are specified according to the tensor strnctnre of the corresponding operators. To distin- 
gnish the nncleon momentnm independent matrix elements (|3^ ) from the momentnm dependent 
ones, (^) and (|3^), the latter will be labeled by an extra index p. AT reflects the isospin se¬ 
lection rule in the nuclear transitions. If the isospin of both the initial and final states equals to 
zero then AT = 0, and if the isospin of one state is T=0 and of the other is T=1 then AT = 1. 
Only when both isospins are equal to 1 we have both values; AT = 0 and 1. 

In this paper we will discuss: 


— the ®Li and ®IIe observables in their ground state, 


— the transition J’^T = 1’''0 —O’*"! and vice versa generated by various projectiles. 


— the 7 -transition to the J'^T =2~^1 level and 


— the transition J'^T = l+O —>• J'^T = 3’''0 . 


Matrix elements involved in the calculations of the corresponding observables are the follow¬ 
ing: 

1. The ground state of ®Li (AT = 0): 

— i) charge and body (matter) rms radii - [ 000 ]o; 

— ii) quadrupole moment - [ 022 ]o; 

— iii) longitudinal form factor - [ 000 ]o and [ 022 ]q; 

— iv) magnetic moment - [101]o and [lll;p]o; 
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— v) magnetic elastic form factor - [101]o, [121]o and [lll;p]o; 

— vi) the elastic pion scattering - [000]o, [022]o, [101]o and [121]o. 

2. The ground state of ®He (AT = 0 and 1): 

— i) rms charge radius - [000]o and [000] i; 

— ii) rms body radius - [000]o- 

3. The transition from the ground state of ®Li to the J'"T = O’*"! level in ®Li or in ®He and 
vice versa: 


— i) beta decay - [101]i; 

— ii) muon capture - [101]i, [121]i, [011;p]i and [lll;p]i; 

— hi) form factor - [101]i, [121]i and [011;p]i; 

— iv) positive pion photoproduction - [101] i and [121] i. 

4. The gamma transition to the J'"T = 2’''1 level of ®Li: 


— i) [101]i and [lll;p]i. 

5. The transition from the ground state of ®Li to the J'^T = 3’''0 level in ®Li: 


— i) gamma transition and electron scattering - [022]o, [044]o and [112; p]o; 


— Ill 


pion scattering - [022]o, [044]o, [12u]o and [14m]q. 


Among the above listed matrix elements only [000]o, [022]o, [044]o and [lwu;p]o are of pure 
spatial type. For overlapping clusters, as has been shown in Ref. ordinarily the magnitude 
of the exchange effects for a pure spatial operator is associated with the quantity 6 = 1 — Q^, 
where is given by eq. (|2D|). The MDM functions [|^], which are used throughout this paper, 
were constructed in a way to minimize the admixture of components with Pauli-forbidden Young- 
schemes [f] = [51] and [6] (see Ref. [^). For this reason the values of are rather close to 
1 and therefore the exchange effects are usually not large. That is why such observables as ®Li 
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charge and body radius, quadrupole moment and ®He body radius should have very close values 
in all versions of the model. The direct calculations have demonstrated (see Table 1) that the 
renormalization effects in this group of matrix elements is indeed small. 

Matrix elements [lOljo and [101]i, in the long-wave approximation where f{qrj) = Const, 
are of pure spin and spin-isospin types, respectively. The kernel renormalization restores exactly 
their values, altered by the C-renormalization. At larger momentum transfer, as follows from 
Table 1, they depend on the renormalization procedure. Now, let us turn to the discussion of 
the renormalization effects on various observables of the six-nucleon nuclei. 

B. The static properties of the six-nucleon systems 

In Table 2 the results of calculation of the static properties of the ®Li (columns 2-4) and ®He 
(columns 5 and 6) are presented. 


1. Ground state of^Li 


Only one matrix element contributes to the rms radius and the quadrupole moment of ®Li. 
Sensitivity of these observables to the version of the renormalization procedure have been dis¬ 
cussed above. Here one can only add that the calculated rms radius is in agreement with the 
experimental data whereas the quadrupole moment is far from the measured one. The theoret¬ 
ical prediction for the value of the quadrupole moment is a special problem and it is beyond 
the scope of this paper. It seems that one should take into account the D-component of the 
a-particle wave function to come close to the experimental value of quadrupole moment of ®Li. 
It is worth mentioning that in all large-scale calculations of this moment (see, for example, the 
discussion in Ref. [^] and references therein) its predicted value is very close to that obtained 
in this paper. 

One can see from Table 2 that both the exchange and renormalization effects for all static 
observables of six-nucleon system are not large. 
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2. Gamma transition to the J'^T = 2^1 level in ^Li 


The exchange effects in the static characteristics can increase if the leading matrix element 
of the corresponding operator is suppressed due to nuclear structure effects. Such a situation 
occurs in ®Li for the Ml 7 -transition from the ground state to the J'^T = 2~^1 level. Here the 
magnetization current matrix element [ 101 ] 1 is strongly suppressed due to the structure of the 
corresponding wave functions (the corresponding operator O is unable to change the total angular 
momentum L by two units) and the main contribution comes from the convective current, [lll;p]. 
So, this transition belongs to the scissor mode. The prediction of the MDM and of two versions 
of the AMDM differ noticeably from each other: 

(eV) = 0.29 (MDM), 0.38 (AMDMc) and 0.34 (AMDMk). The experimental value is 
(0.27 ±0.05) eV. 

The result obtained within the framework of the MDM is closer to that of the AMDMk, 
than to the AMDMc one. The AMDMk result deviates a little from the measured value. But 
it is expected that in this hindered transition, the mesonic exchange currents could play a very 
important role. So it will be interesting to continue to investigate this problem further. Notice 
that for many other observables the AMDMk results are also close to those of the MDM. That 
is why in these cases the MDM has been very successful in describing these observables. 


3. Ground state of^He 

The neutron halo nucleus ®He is one of the interesting peculiarities of the six nucleon system. 
The outer neutron pair, forming this halo, extends to large distances and is localized in the so 
called quasiasymptotical region. The charge distribution in ®He is associated with two types of 
matrix elements - [ 000]0 and [ 000 ] 1 , contrary to the matter distribution which is associated only 
with the isoscalar part. So, we can compare the results obtained for two similar operators in the 
coordinate space which, however, have different isospin structure. This comparison is done in 
Table 2. 

When the kernel renormalization is adopted within the AMDM, the rms radius of ®He is 


changed by less than 1% when compared with the MDM case. Let’s dehne, following Ref. [33], 
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the size of the halo in ®He as the difference between the proton and neutron radii r^aio =< '^n 
— < Tp It changes by about 7% when one goes from the MDM to the AMDMk despite its 
quasiasymtotic nature. The exchange and renormalization effects for the halo in ®He are not well 
pronounced. As a rule one should keep in mind this effect when calculating the radius of the halo 
in heavier nuclei, where the antisymmetrization between nucleons located in different clusters is 
neglected. More complicated intrinsic structure of the constituent cluster in such systems, ®Li 
cluster in ^^Li as an example, would increase the exchange effects and, as a result, would increase 
the calculated halo radius. 

The antisymmetrization procedure with renormalization to a constant distorts the asymptotic 
behavior of the wave function. As a result a nonphysical exchange effect at large distances arises 
in this case. This is the reason why the body rms radius of ^He changes when going from 
the MDM to the AMDMc. To avoid this unphysical effect one should deal with the kernel 
renormalization from the beginning. 


4 . Beta decay of^He to ^Li 


The ®He /3-decay rate is determined by the spin-isospin Gamow-Teller operator dr: 


ft = 


2/to^ 


+ 0 + 


(38) 


(WsvO'iAifr 

We use the following values for the constants in eq. (P^: 2/fo+^o+ = 6144 s and {qa/ 9v) = 
— 1.259. The results of the ft calculation are the following: ft = 796 s (MDM and AMDMk) 
and 772 s (AMDMc). The experimental value of ft is 813 s. Keeping in mind that there should 
be some meson exchange current contribution, it seems that the calculated value of the /3-decay 
matrix element is somewhat overestimated. The AMDMk result is exactly equals to the MDM 
one, as it should be. 


C. Electromagnetic form factors 
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1. Longitudinal form factors 


The results of the calculations of the longitudinal elastic form factor of ®Li are given in Fig.2. 
Starting from this figure and in all subsequent ones we present the results in the following way. 
By the solid line we show the results obtained within the framework of the AMDMk, by the 
dashed line — within the AMDMc and by the dotted line — within the MDM. The full elastic 
form factor consists of the sum of the contributions of the CO and C2 terms. However the C2 
form factor is visible only in the region of the CO minima which is located at about 3 fm~^. That 
is why we are not showing them separately. All longitudinal form factors for transitions between 
nuclear levels with isospin T = 0 are associated with pure spatial operators. At low momentum 
transfer the elastic charge form factor of ®Li is determined by its rms charge radius which is 
practically not affected by the nucleon exchange. More importantly, the exchange effects remain 
very small up to a high momentum transfer. The effects of the kernel renormalization are too 
small to be seen in the Fig.2. So, one can safely use the MDM for this elastic longitudinal form 
factor calculations. 


As it had been discussed earlier in Ref. and confirmed by the present calculation, in order 
to reproduce the ®Li elastic form factor it is very important to use the correlated a-cluster wave 
function which reproduces the a form factor. Otherwise, the calculated high momentum part of 
the ®Li elastic form factor will be much lower than the measured one. 

The antisymmetrization results in the mixing of the components of the three-body functions 
which have different partial angular momenta I and A but the same total L. In the ®Li ground state 
wave function this situation is not well pronounced, because the dominant weight is associated 
with the I = X = L = 0 component. In the wave function of the J'^T = 3+0 level in ®Li there are 
already two such components: Z = 2, A = 0 and Z = 0, A = 2 with the total angular momentum 
L = 2 and spin S = 1. The antisymmetrization with renormalization results in their additional 
mixing, which differs essentially from the MDM mixing. For the same reason there is noticeable 
difference in the wave functions of this level in the two versions of the AMDM - with and without 
kernel renormalization. Thus the modification of this inelastic longitudinal form factor due to 
the antisymmetrization and renormalization should be rather large. 
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Fig.3 displays the E2 )—the square root of the reduced quadrupole transition strength 

to this level at low momentum transfer region. Just this multipole dominates in the form factor for 
this transition. The kernel renormalization effects are visible already at zero momentum transfer 
due to this additional mixing. At photon point B{E2,1~^0 —>■ 3’''0) in equals to 20.05, 

15.25 and 20.71 for the MDM, AMDMc and AMDMk models, respectively. The experimental 


value is 21.8 ± 4.8 |^. So, more precise measurement of this observable seems to be needed to 
test the theoretical models. Fig.4 displays the form factor for this transition in a broad region of 
momentum transfer. In logarithmic scale, it is not easy to see the difference between the three 
versions of the model. In the region of the hrst maximum this difference is about the same as in 
the region of very low momentum transfer. The magnitude of the calculated form factor is lower 
than the experimental one. The reason for this is not clear enough. At high momentum transfer, 
for example, it is possible that noncentral components of the a-particle wave function, which are 
not taken into account at all, start to play a role. So, this discrepancy should be a subject of 
further studies. 


2. Magnetic form factors 


The earlier MDM calculations of the ®Li magnetic form factors in Ref. |^] have demonstrated 
quite clearly that despite of a good agreement with experimental data in the region of the hrst 
maximum, the theoretical result for the region of the second maximum lies much below the 


experimental data for any choice of NN and a-N interaction [40,41]. At the same time, the 


studies carried out in Ref. [^] showed qualitatively that the antisymmetrization is of great 
importance when one calculates the ®Li transversal form factors. Both the elastic and the inelastic 
transversal form factors to the J'^T = 0^1 level are governed by the Ml multipole only. The 
antisymmetrization of the wave function leads to a strong increase of the form factor within the 
region of the second maximum, thereby improving the agreement with the experimental data 
substantially. This is demonstrated in Fig.5 and Fig.6. 

At low momentum transfer the Ml operator is proportional either to the operator of the 
magnetic moment (in the case of the elastic scattering) or to the GT-operator of /J-decay and has 
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mainly a spin-isospin component. Here we use the fact that the convective current contribution 
is small. That is why the exchange effects are very small in this region of momentum transfer. 
However, already at the second maximum their role becomes crucial because the spatial part of 
the Ml operator starts to play its decisive role. 

It is necessary to stress that in the AMDMk version the position of the minimum is reproduced 
very precisely. The other interesting effect of kernel renormalization consists in stabilization of 
the AMDM results for intermediate momentum transfer (g < 2.5 — 3.0/m“^): the magnitude 
of the form factors in the region of the second maximum becomes almost independent on the 
details of the a-cluster wave function. It means that in this region the exchange and short-range 
correlation effects are well separated from each other. This result is very important because our 
a-cluster wave function has a pure phenomenological origin. 

At larger momentum transfer an interference between two types of nucleon correlations (ex¬ 
change and short-range) starts to play its decisive role. The consistent treatment of both effects 
results in appearance of a third maximum in the Ml form factors at about ~ 20/m“^. Again 
as it happened with the second maximum, the AMDMk version provides more stability of the 
third maximum in the magnetic form factor as compared to the AMDMc version. 

Elastic Ml form factor has an isoscalar nature. Therefore, the exchange currents in this 
case are not very strong and should not alter considerably the results for small and moderate 
momentum transfer (see discussion in Ref. f^). The inelastic Ml form factor has an isovector 
origin and it is somewhat lower than the experimental data, already at lower momentum transfer. 
One can assume that this is a real indication for the mesonic exchange currents in this transition. 
At higher momentum transfer their contribution should increase. 

It seems to us that it is very promising to continue the experimental study of ®Li form factors 
and to cover the high momentum transfer region. This can be an interesting topic, for example, 
in CEBAF. 
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D. Gamow-Teller transition in muon capture 


The Gamow-Teller transitions in muon capture on nuclei are interesting from the following 
points of view: 

- how well the nuclear models are able to predict their strength at intermediate momentum 
transfer, 

-what is the magnitude of the induced pseudoscalar coupling gp m. nuclei, 

-how large are the mesonic exchange currents at this momentum transfer. 

Up to now these problems have been analyzed in great details only in the few-body systems. 
®Li gives the other opportunity to study the above discussed problems in a nucleus with twice as 
many nucleons as ^He. Here, all these effects can be enhanced. 

The momentum transfer in muon capture on ®Li with formation of ®He in its ground state 
is equal to |^ = E^/c = 100.7 MeV/c. This value is in the region of the hrst maximum in the 
inelastic electron scattering form factor (see Fig.6) to the T = 0^ 1 level of ®Li. Four matrix 
elements form two transition amplitudes, Ti and T 2 , to the ground state of ®He: 

ri = 17dM{i_Ari} (39) 

and 


P = (2/9)'-'"(G^ - Gp)[101|{l + AP], 


(40) 


Here we have isolated the dominant matrix element [101] from the rest. ATi and AT 2 are built 
up from the three others, [121], [11 Ip] and [01 Ip]: 


and 




AT, = 


(41) 


(42) 


[ 101 ] Ga — Gp iV7;v[101] 

Ga and Gp are some combinations of the weak interaction coupling constants. They are given 
in Appendix A together with the expression for the capture rate. 
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Two observables are of experimental interest in this partial muon capture transition - the 
capture rate A, when the hyperhne states are populated according to their statistical weights, 
and the ratio of the capture rates from the upper, A+ , and the lower, A_ , hyperhne states of 
the ®Li mesic atom. 

The capture rate is proportional to the sum of Tf and the ratio A_|_/A_ is given by the 
following expression: 

^ {3Ga - Gp}^ ^ 

where 


A2 


3{Gp — Ga) L 3 
2G^ I ~ ^ 



(44) 


and 


ATi + AT2 
1 +AT 2 


(45) 


It is through X that the nuclear structure is manifested in the ratio of the capture rates A+/A_. 
In the limit ATi and AT 2 —> 0, X and A 2 are equal to zero and this ratio is given by the coupling 
constants only showing a strong dependence on gp. That is why this ratio is a subject of detailed 
analysis in many nuclei (see Ref. ji^). 

ATi is independent of gp. After hxing the weak interaction constants according to Appendix 
A one arrives to the values of Ti and ATi listed in Table 3. As follows from Table 3, the 
combination of matrix elements which forms ATi is sensitive to the nuclear model. But again, 
the result in the AMDMk version is closer to the MDM one. The difference is about 10%. The 
full amplitude Ti is less sensitive to the model because ATi itself is about half of a per cent. The 
combination of the matrix elements which forms AT 2 is not very sensitive to the model as well. 
But its relative weight in the full amplitude T 2 is about 15 %. 

In both cases it is important to take into account the velocity dependent matrix elements 
(the corresponding operator is proportional to the momentum of nucleon inside the nucleus) . If 
one neglects them (see last column of Table 3), as it is done in some papers, then both terms, 
ATi and AT 2 , will be changed drastically. 
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The ratio of the hyperhne capture rates depends on the nuclear model through X sensitivity of 
the amplitudes and their combinations to the nuclear model we give in Table 3 their magnitudes 
for hxed value of the induced pseudoscalar coupling constant {gp/Q a = 7). Finally, for the same 


value of the coupling constant we give the capture rate and the ratio A+/A_. 

Fig.7 demonstrates the calculated capture rate as a function of Qp/qa together with experi¬ 
mental data [^]. They are far from each other. Even if one increases arbitrarily the theoretical 
capture rate by about 10% , to take into account in this way the meson exchange current contribu¬ 
tion, the situation improves only slightly. But nevertheless, the value of the induced pseudoscalar 
constant which is needed to reproduce the experimental data is lower than that which follows 
from the Goldberger-Treiman relation. It seems very desirable to remeasure this capture rate 
with a better accuracy to insure that situation is not due to the existing experimental data. 

Fig.8 demonstrates the predicted values for the ratio of the capture rates as a function of 
Qp/qa for the three versions of the model. The absolute value of this ratio appears to be the 
same within about 3%. It means that the uncertainty from the nuclear structure side is minor. 
As the next step, one should incorporate into the calculations the meson exchange currents, 
which can modify the result obtained for the ratio of the capture rates. 


E. Partial transitions in pion scattering and pion photoprodnction off ®Li 

After successfully employing the developed model for description of the ®Li traditional ob¬ 
servables let us turn our attention to the: 

— pion elastic and inelastic scattering with the excitation of the J'^T = 3’'“0 level, 

— pion photoprodnction with the formation of ®He in its ground state (J’^T = 0+1). 

These reactions are selected for the following reasons. Pion scattering on polarized target has 
already been carried out and the vector analyzing power iTn for the ground state and the 
J'^T = 3+0 level transitions was measured. Nuclear models can be tested in a more thorough 
way via the polarization observables than via the differential cross sections, provided that the 
reaction mechanisms are properly taken into account. So, pion scattering can be considered as 


29 



complementary to electron scattering in probing details of the ®Li structure. One can invert 
the formulation of the problem — if the nuclear structure is well determined, then the other 
ingredients of the theory can be checked. 

As to the pion partial photoproduction off nuclei it is expected that very soon the data with 
polarized photons will be available as well. It will be data on the beam asymmetry: 


m = 


a^{9)-a^9) 


(46) 


(t-^{9) + 0 - 11 ( 6 ') ’ 

where 9 is the polar angle of the outgoing pion, is the differential cross section for the photon 
polarization along the x axes, and — along the y axes. The pion momentum q is directed 
along the x axes, the photon momentum k along the z axes and the vector [k x q\ along the y 
axes. We will demonstrate that, depending on the momentum transfer (pion emission angle), 
there are regions sensitive either to the nuclear model or to the details of the reaction mechanism. 
At present there are polarization measurements in pion photoproduction when nucleon knocked 
out from nucleus . Pion scattering and pion photoproduction on light nuclei were described very 
successfully within the framework of DWIA in momentum space (see for example Ref. P5| , ^ ). 
However, in some cases (including elastic and inelastic to the 3’'“0 level scattering in ®Li) one has 
to use the coupled channel method (see for example Ref. [P7| ). 

In this paper we will continue to exploit the coupled channel method for pion scattering to 


the ground state, the J'^T = 1’''0 and the 3’''0 levels, in line with Ref. E^. Pion photoproduction 


will be analyzed within the framework of the DWIA, in line with Ref. |^6[|. Some basic formulas 
for both processes are given in Appendix B. 

Among all low lying states of ®Li the elastic and the inelastic to the 3’''0 cross sections are 
measured for several energies. Cross sections to these levels are large due to the isoscalar nature 
of the nuclear transitions. More or less complete data, which include the polarization observables, 
exist for a few energies of the incoming pions. One of such energy is = 134 MeV. That is why 
we will concentrate all our attention mainly on the data at this energy. 
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1. Elastic and inelastic to the J'" = level pion scattering 


Elastic and inelastic to the 3'''0 level pion scattering off ®Li are governed by the spin indepen¬ 
dent part of the elementary ttN scattering amplitnde. The spin-dependent part of pion-nncleon 
amplitnde is mnch smaller and does not show np in the differential cross section. In this section 
we compare onr resnlts with the experimental data of Ref. . 

Fig.9 displays the elastic cross section for 134 MeV pions. All three versions of the model 
give similar results. The situation is the same as in the case of electron scattering where the 
calculated longitudinal form factors appeared very close to each other. The point is that in pion 
scattering at this energy the maximal momentum transfer does not exceed 2.5 fm~^ where these 
three versions of the model predict very close values for the matrix element — see Table 1. 

In the case of the 3^0 level the AMDMc version at low momentum transfer region (small 
angles) gives the lowest result just as in the electron scattering case (see Fig.10). The MDM and 
AMDMk versions predict a higher value which is more close to the experimental data. Predictions 
of all three versions of the model in the region of the minimum differ from each other significantly. 

In both cases one can say that agreement between theory and experiment is quite good and 
some deviation of the theory from experiment can be due to some shortcomings associated with 
the reaction part. For example in this calculation the pion absorption is ignored. Of course this 
process can not play a decisive role, but it can have some effects. 

Contrary to the differential cross section the vector polarization Tn is caused by interference 
between the spin-dependent and the spin-independent parts of the pion-nncleon amplitude. To 
have some feeling of Tn it is instructive to use the plane wave approximation, following Ref. 


In the elastic longitudinal electron scattering form factor one could safely neglect the L=2 matrix 
elements. The same is true for the pion scattering. Then in the plane wave approximation one 
gets 


iTn{da/dn) ~ [000][101]Jm[/^Ar(s = 0)/:^(s = 1)]. (47) 

Here fnN{s = 0) and firNis = 1) are the spin-independent and spin-dependent parts of the 
elementary iiN scattering amplitudes, respectively. From this expression it follows that iTn 
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changes its sign where the Ml form factor does. The largest value of the vector analyzing 
power corresponds to the minimum in the differential cross section of pion scattering. The 
Im[f.,rNis = 0)/*^(s = 1)] changes its sign at 6 larger than 120° for > 120 MeV. Thus the 
vector analyzing power is a complicated function of the scattering angle. 

The vector analyzing power iTn for the transition to the 3’''0 level in the plane wave approx¬ 
imation looks very similar: 

tTn{da/dn) ~ [022]{[122] + (32/49)[123]}/m[/^^(s = 0)/;^(s = 1)]. (48) 

The results of comparison between experimental data and the three versions of the model is 
given in Fig. 11 for the elastic scattering and in Fig. 12 for the transition to J'^T = 3’''0 level. 

Contrary to the differential cross section the vector analyzing power is more sensitive to the 
nuclear model. All three versions of the model used in the calculations give different results. The 
AMDMk version explains the shape of the vector polarization as a function of scattering angle. 
There is some discrepancy in the magnitude of the vector analyzing power. This discrepancy can 
be caused by the approximate treatment of the pion-nuclear dynamic. But what emerges from 
the calculations is that the nuclear model is already reliable enough so one can begin to improve 
upon the dynamic. Simultaneously it seems to be very useful to have more precise experimental 
data on this quantity. 


F. Pion photoproduction to the ®He ground state. 


Pion photoproduction at threshold is dominated by the Kroll-Ruderman spin-flip term. So 
in the pion photoproduction on ®Li with formation of ®He in its ground state the same matrix 
elements governing the transition as in muon capture and in electron scattering. So at threshold 
one expects a successful description of experimental data. In Fig. 13 the calculated results within 
the framework of the AMDMk are compared to experimental data of Ref. P9| , |^ at 200 MeV 
energy of the incoming photons. At this energy the momentum transfer is covered from 0.5 fm~^ 
at 25° to 1.5 fm~^ at 135°. At this end point pion photoproduction is governed by the form 
factor of electron scattering in its minima, where it is very difficult to expect a precise theoretical 
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description of the momentum dependence. That is why at backward angles there is a serious 
disagreement between theory and experiment. As to the forward angles one sees a nice agreement 
with experimental data. 

At higher energies of incoming photons, say at 320 MeV, contribution of A-isobar becomes 
very important and it brings the spin-independent part into the amplitude. At this energy the 
difference in the differential cross section between the three versions of the model is very small 
at forward angles which correspond to low momentum transfer. Fig. 14 demonstrates this in a 
clear way. Predictions become different from each other starting from the angle 6 at about 60°. 

The beam asymmetry at forward angles is also weekly sensitive to the model — see Fig. 15. 
At these angles a strong sensitivity of the beam asymmetry to the A-isobar property inside a 
nucleus is found. In Fig. 16 the beam asymmetry is plotted as a function of the A-isobar mass, 
Ma, in nuclei. A 5% deviation of this mass from its free value changes the beam asymmetry 
noticeably. This happens in the region of angles from 20 to 40° where the cross section is not 
small. So, this result can be considered as a motivation for measuring the beam asymmetry in 
this transition. Variation of the A-isobar mass affects the differential cross section mainly at 
backward angles (Fig. 17). However at these angles there is a strong sensitivity to the nuclear 
model as well. 


IV. CONCLUSION 

A modihed version of the Multicluster Dynamic Model was proposed to construct completely 
antisymmetrized wave functions of multicluster systems. In the model the intrinsic wave functions 
of the constituent clusters were inserted explicitly, and as a result the model became closer to the 
MCM, where the motion of all nucleons is taken into account on equal ground. However, contrary 
to the MCM, the cluster-cluster potentials used are those constructed on the basis of phase shifts 
analysis of the free cluster scattering and thus the on-shell properties of these potentials are 
taken into account in the most precise way. Such a modihcation, nevertheless, retains the easy 
application of the MDM and at the same time allows us to incorporate into the calculations the 
realistic wave functions which describe the motion of the nucleons inside the clusters. 
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A group-theoretical analysis of matrix elements for arbitrary multicluster systems was de¬ 
veloped. On its basis all operators of single-particle origin have been classihed according to 
their response to the antisymmetrization procedure. A group of operators is found whose ma¬ 
trix elements are strictly conserved when going from the MDM to the AMDMk, but not to the 
AMDMc. This analysis allows to predict, before the numerical calculations are carried out, how 
essential will the exchange effects be for the particular nuclear characteristic. Simultaneously, 
it was shown that the kernel renormalization removes the spurious effects which had appeared 
when the renormalization to constant was adopted. 

Detailed numerical calculations of various observables in the six nucleon systems were per¬ 
formed. These nnclei are particularly singled out in nnclear physics becanse they already exhibit 
qnite properly the typical properties of all light nuclei, but at the same time the number of 
nucleons is not too large to prevent an accnrate and consistent calculations starting from the 
realistic interactions. 

Within the framework of the AMDMk a better description of many six-nncleon observables 
has been achieved. The electron scattering form factors are those obtained with high precision. 
One of the interesting effects of the kernel renormalization consists in the stabilization of the 
AMDM results for intermediate momentum transfer (g < 2.5 — 3.0 the magnitude of the 
magnetic form factors in the region of the second maximnm became almost independent on the 
details of the a-cluster wave fnnction. It means that in this region the exchange and short-range 
correlation effects are well separated from each other. This result is very important for the 
accnracy of the model because our a-cluster wave function has a pure phenomenological origin. 

At larger momentum transfer an interference between two types of nucleon correlations (ex¬ 
change and short-range) starts to play its decisive role. The consistent treatment of both effects 
results in the appearance of a third maximum in the Ml form factors at ~ 20fm~^. And 
again the AMDM^- version provides some stabilization of this maximnm. 

The model developed was applied to mnon capture, pion scattering and pion photoproduction. 
It succeeded in providing a good description of the differential cross sections in pion scattering. 
What is more important, it allowed us to reproduce at least the shape of the vector polarization 
which is more sensitive to the details of nuclear structure. Thus we have demonstrated that the 
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AMDMi^ allows us to accurately describe tiny details of the nuclear structure which, nevertheless, 
are essential for the interpretation of the polarization phenomena in nuclear reactions. 

In the case of pion photoproduction, a region of angles was found where the beam asymmetry 
depends only slightly on the nuclear model and very strongly on the mass of the A-isobar inside 
nuclei. Experimental data in this region of angles, if they appear, can provide fresh insight into 
this problem. 

The calculations performed have demonstrated that the results of the group-theoretical analy¬ 
sis are valid in a more realistic case as well, when the cluster wave functions are more complicated 
as compared to the SU(4) and SU(3) scalars. 

There is one more aspect of this paper which can be formulated in the following statement: 
after introducing the kernel renormalization the model became more universal and can be applied 
now to more heavy multicluster systems, contrary to the MDM. 

The group-theoretical analysis carried out in this paper offers promise as a processing tool 
and can be used in many other situations. Here we would like to demonstrate two more areas 
of its application. For a certain (and rather wide) type of light systems the many nucleon wave 
functions constructed within the framework of the oscillator shell model can be rewritten in bi- 
or multicluster [^] representation. In other words such a wave function can be represented 
by eq.(]T2|) where the a-particle, etc., and nucleons act as clusters. So the group theoretical 
analysis may be applied to such systems as well. From this analysis one is able to get early 
insight on the role of the nucleon exchange in the discussed systems and to estimate its effect 
arising from the antisymmetrization between nucleons in different clusters. We will not discuss 
all these issues because they are beyond the scope of this paper. 

There is another interesting held of application for the developed scheme of the group- 
theoretical analysis. The case in point is the problem known as quarks inside the nuclei. If 
one starts from the three-quark representation of the individual nucleons in nuclei (and does not 
take into account contribution of ”sea”), then the question whether one should take into account 
the quark degrees of freedom of nuclei, or it would be enough to consider only the nucleons as 
constituents of nuclei, depends substantially on the magnitude of the quark exchange effects. 

Indeed, as we have seen, in the multicluster system any process without the breaking down 


m 
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of the constituent clusters can be described within the framework of the approaches where the 
clusters themselves are treated as an elementary particle. Such an approach is valid when the 
exchange effects between clusters are small. 


To estimate the quark exchange effects in nucleus it is again convenient to use the methods 
of group-theoretical analysis. Let us discuss very qualitatively two conclusions which follow from 
such an analysis. 

1. The exchange effects should be larger for observables which are described by the full tensor 
operator. However even in this case the quark exchange effects in nuclei are usually rather small 
because of high internucleon distances. That is why the quark degrees of freedom are seen in 
such processes only at very high momentum transfer. 

2. Simultaneously from our analysis (by reduction ab absurdum) it follows that the quark 
effects in nuclei should be more pronounced in processes where the break down of the quark 
clusters or their excitation take place. Such processes are enhanced due to the exchange effects 
0. One interesting example is already considered in Ref. |5^, where the nucleon resonance as 
a spectator participates in the reaction. Though the weight of the components with the nucleon 
resonance in the nuclear ground state wave function is very small, observation and investigation 
of the corresponding reaction is important due to more transparent manifestation of the quark 
degrees of freedom in nuclei in this case. 

Thus, we have shown in this paper that the group-theoretical analysis can open a new way 
for general investigation of exchange effects in various nuclear and subnuclear processes. Some 


specihc examples have been discussed, for example, in Ref. |54 


ACKNOWLEDGMENTS 

This work is supported in part by the National Science Foundation under Grant ^ PHY- 
9413872 and Russian Foundation of Basic Research Grants 96-02-18691. 

R.A.Eramzhyan would like to express his thanks to the Gyclotron Institute at Texas A&M 
University for the kind hospitality. 

The authors would like to thank A.A.Ghumbalov and S.S.Kamalov for cooperation in the ini- 


36 











tial stage of this work concerning the pion physics, S.Shlomo for careful reading our manuscript 
and many valuable suggestions and V.I.Kukulin and V.N.Tolstoy for the numerous fruitful dis¬ 
cussions. 


APPENDIX A: GENERAL EXPRESSIONS FOR MUON CAPTURE RATES 


The muon capture rate is expressed via the amplitudes given by eqs.(^) and in the 
following way: 

1 {2Jf + 1) 


A = STi'^a'^{mfj,/mejln2ri{Ey)Z'^Rz 


2/to+^o+ (2^1 + 1) 


{t2 + t|}. 


(Al) 


Here 


ri{Ey) = {Ey/m^f 1 


E,, 


m,. 


1 3 


+ Mi j l + m^/Mi_ 


(A2) 


Rz = 0.92 is the dimensionless quantity and is equal to the ratio between the averaged muon 
wave function squared over the nuclear volume and its value for the point nucleus at the origin, 
Z is the charge of ®Li and Mj is its mass. is the nucleon mass. 

After eliminating the main matrix element [101] from the amplitudes Ti and T 2 one comes to 
the following expression for the sum of Tf and 

n2 I rr2 2Gq_2^[101]^ 


Tf = 


-{1 + Ai} 


(A3) 


where 


Gl_j. = G\ + Gl/3 - 2 GaGp/3 


is the square of the Gamow-Teller constant in muon capture, 

(Ga - GpY 


A 


{AT2-2ATi} 


+ 


3 Gq_j. 


{ATI + 2AT2} 


and 


(A4) 


(AS) 


GA{q ) = 9A{q ) -4:.706gv{q ) 


2\ 


2M]\[ 


(A6) 
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Gp{q^) = {gp{q^) - gA{q^) - ^Gmgy{q^)} 


Ev 

2M]y 


(A7) 


where gA is the axial-vector coupling constant, 


gA{q^)Jgv{^) = -1.239, 


(AS) 


and 


gA{Q)/gv{Q) = -1.259. 


(A9) 


gv is the vector coupling constant 


9v(9")/9v(0) = 0.972, 


(AlO) 


and gp is the induced pseudoscalar coupling constant, which is less known compared to the two 


others. It is hxed by the Goldhaber-Treiman relation to be 


gp{q^)/gA{Q) = 6.78. 


(All) 


Ordinary calculation results are given versus this coupling constant. 


APPENDIX B: GENERAL FORMALISM FOR PION SCATTERING AND PION 

PHOTOPRODUCTION 

The starting point for pion scattering and pion photoproduction calculations is the con¬ 
struction of the multiple scattering matrix T(E). It is obtained as a solution of the Lippmann- 
Schwinger equation 

T{E) = V{E) + V{E)PG{E)T{E). (B1) 

Here G(E) is the pion-nucleus Green’s function, V(E) is a potential matrix, and P is a projection 
operator. Projection to the nuclear ground state corresponds to application of the optical model, 
whereas projection into a group of nuclear states corresponds to the application of the coupled 
channel method. 
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The potential matrix is usually divided into a microscopic first order term Vi{E) and a 
phenomenological second order term V 2 {E) : V(E) = Vi(E) + V 2 (E). The potential matrix Vi{E) 
contains the full spin and isospin dependence of the pion nucleon amplitudes via the impulse 
approximation and receives contributions from the nuclear matrix elements. The second order 
term is associated with true pion absorption and higher order processes. As follows from Ref. 

4^ this term is weak in ®Li. To simplify the calculations we neglect this term. 

The pion photoproduction off nuclei is treated by the DWIA method. The pion photopro¬ 


duction matrix is obtained from the equation similar to that in eq.f ^ll) : 


T^^ = U^^ + Tl^,G{E)U., 


TT'Y. 


(B2) 


Here is the plane wave photoproduction amplitude off nuclei, 

E{q) = ET^{q) + E^iq) is the total energy of the pion-nuclear system. The auxiliary matrix 
T' is related to the pion-nuclear scattering T-matrix discussed above by the following way: 

T'(g,go) = p-l)M]T(g,g„). (B3) 


In the Impulse Approximation the matrix elements of U.,^^ are expressed in terms of the elementary 
(7,7r) amplitude. In our calculations the Blomqwist-Laget amplitude for pion photoproduction 
(see Ref. |^) is used. 
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TABLES 


TABLE 1. Renormalization effects for some matrix elements in electron scattering on ^Li. A — for 
elastic scattering, B — for the transition to the J'"T = 3^0 level, C — for transition to the J^T = 0+1 
level 



MDM 

AMDMc 

AMDMk 

q = 0.5/m ^ 

q = 2.5/m ^ 

q = 0.5/m ^ 

q = 2.5/m ^ 

q = 0.5/m ^ 

q = 2.5/m ^ 

A 

O 

O 

o 

o 

.4669 

6.290E-3 

.4723 

7.981E-3 

.4657 

5.261E-3 

B 

[022]o 

3.825E~2 

1.621E-2 

3.375E-2 

1.752E-2 

3.888E-2 

1.656E-2 


[101]i 

-.2004 

4.235E-3 

-.2081 

1.081E-2 

-.1954 

8.840E-3 


[121]i 

-3.839E-3 

-5.728E-4 

-4.016E-3 

-8.531E-4 

-4.152E-2 

-7.007E-4 

C 

[lllp]i/Mjv 

-2.199E-3 

-4.629E-4 

-2.883E-3 

-6.968E-4 

-2.684E-3 

-6.043E-4 


[011p]i/MAr 

7.303E~3 

-1.097E~3 

7.731E-3 

-1.939E-3 

7.728E-3 

-1.546E-3 


TABLE 11. Static properties of six-nucleon nuclei 


Nucleus 

®Li 

h/ho 

®Li 

Q, fm^ 

®Li 

®He 

rhaioJm 

®He 

MDM 

2.55 

0.829 

0.49 

2.43 

0.74 

AMDMc 

2.48 

0.838 

0.49 

2.33 


AMDMk 

2.55 

0.829 

0.51 

2.44 

0.80 

Exper. 

2.56 

0.822 

-0.082 

2.33±0.04 

0.87±0.06 

32] 
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TABLE III. Amplitudes of muon capture in ^Li, their ratio and capture rates for three versions of 
the model. Numbers with * are given for results obtained without the velocity terms. 



MDM 

AMDMc 

AMDMk 

ATi 

4.74E-3 

2.47E-3 

3.97E-3 




1.56E-3* 

ATs 

-.136 

-.139 

-.147 




0.031* 

Ti 

0.198 

0.205 

0.193 

T2 

0.0687 

0.0708 

0.0660 

X 

-0.152 

-0.159 

-0.168 

A, s—1 

1259 

1351 

1192 




1225* 

A+/A- 

0.0417 

0.0423 

0.0431 




0.0247* 
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Figure captions 


Fig.l. The Jacoby coordinates of the a-2N system. 

Fig.2. Longitudinal elastic form factor of ®Li. The solid line — calculation within the frame¬ 
work of the AMDMk, dashed line — within the AMDMc and dotted line — within the MDM. 
Experimental data from Ref. ||35|| . 

Fig.3. The square root of the reduced quadrupole transition strength B(C2,q) in units efrri^ 
for transition to the J'^T = 3’''0 level in ®Li. Experimental data from Refs. p6|-|38| The notations 
are as in Fig.2. 

Fig.4. The form factor for transition to the J'^T = 3’''0 level of ®Li. Open triangles are 


experimental data from Ref. solid triangles are data from Ref. [Q. The notations are as in 

Fig.2. 

Fig.5. Elastic magnetic form factor of ®Li. Experimental data from Ref. [^. The notations 
are as in Fig.2. 

Fig.6. Inelastic magnetic form factor for transition to the J^T = O’*"! level of ®Li. Experi¬ 


mental data from Ref. The notations are as in Fig.2. 

Fig.7. Muon capture rate for transition to the ground state of ®He. Experimental data from 
Ref. [^. The notations are as in Fig.2. 

Fig.8. Ratio of the hyperhne muon capture rates to the ground state of ®He. The notations 
are as in Fig.2. 

Fig.9. The differential pion elastic cross section on ®Li at =134 MeV. The second order 
pion-nucleus potential is omitted. The notations are as in Fig.2. 

Fig.10. The differential inelastic pion cross section to the J^T = 3’''0 level at T.^ = 134 MeV. 
The second order pion-nucleus potential is omitted. The notations are as in Fig.2. 

Fig. 11. The vector analyzing power for pion elastic scattering on polarized ®Li at = 134 
MeV. The second order pion-nucleus potential is omitted. The notations are as in Fig.2. 

Fig.12. The vector analysing power for pion inelastic scattering to the J'^T = 3+0 level on 
polarized ®Li at =134 MeV. The second order pion-nucleus potential is omitted. The notations 
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are as in Fig.2. 

Fig.13. Experimental and calculated differential cross section for pion photoproduction on 
®Li at Ty = 200 MeV with formation of ®He in its ground state for three versions of the model. 
The notations are as in Fig.2. Experimental data are from Refs. US- 

Fig.14. Calculated differential cross section for pion photoproduction on ®Li at = 320 MeV 
with formation of ®He in its ground state for three versions of the model. The notations are as 
in Fig.2. 

Fig.15. Calculated beam asymmetry for pion photoproduction on ®Li at = 320 MeV with 
formation of ®He in its ground state for three versions of the model. The notations are as in 

Fig.2. 

Fig.16. Calculated beam asymmetry for pion photoproduction on ®Li at = 320 MeV with 
formation of ®He in its ground state for three values of A — isobar mass in nucleus: solid line — 
the mass is equal to the mass of the free particle (Ma), short dashed line — the mass is larger 
by 5% than the free mass, long dashed line — the mass is smaller by 5% than the free mass. The 
calculation is within the framework of the AMDMk. 

Fig.17. Calculated differential cross section for pion photoproduction on ®Li at = 320 MeV 
with the formation of ®He in its ground state for three values of A — isobar mass in nucleus: 
solid line — the mass is equal to the mass of the free particle (Ma), short dashed line — the 
mass is larger by 5% than the free mass, long dashed line — t he mass is smaller by 5% than the 
free mass. The calculation is within the framework of the AMDMk. 
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